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The influence of gas expulsion and initial mass- segregation 
on the stellar mass-function of globular star clusters 
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ABSTRACT 

Recently iDe Marchi et all (|2007f ) studied a sample of twenty globular clusters and 
found that all clusters with high concentrations have steep stellar mass- functions while 
clusters with low concentration have comparatively shallow mass-functions. No glob- 
ular clusters were found with a flat mass-function and high concentration. This seems 
curious since more concentrated star clusters are believed to be dynamically more 
evolved and should have lost more low-mass stars via evaporation, which would result 
in a shallower mass- function in the low-mass part. 

We show that this effect can be explained by residual-gas expulsion from ini- 
tially mass-segregated star clusters, and is enhanced further through unresolved bina- 
ries. If gas expulsion is the correct mechanism to produce the observed trend in the 
c — a— plane, then observation of these parameters would allow to constrain cluster 
starting conditions such as star formation efficiency and the time-scale of gas expul- 
sion. 

Key words: globular clusters: general - stars: luminosity function, mass function - 
stellar dynamics methods: N-body simulations - binaries: general 



1 INTRODUCTION 

The stellar initial mass function (IMF) is one of the most 
fundamental distribution functions of astrophysics. The 
number of stars that form in one event with stellar masses 
in the interval m,m + dm is dN = £(m)dm, where £(m) 
is the IMF. Star- formation theory predicts systematic vari- 
ations of t he IMF with chang e s of the phy s ical condi- 
tions (e.g. iMurrav fc Linl 11999: iLarsonl Il998l: lElmegreenl 



12004 iTillev fe Pudrit dl2005l : iBonnell et al 



20071 ). Neverthe- 



less, the IMF has been found to be essentially invariant 
l|Motte et alJll998l . l200ll : lKroupall200ll . liool ), which is still 

a ch allenge for theoreticians. 

IDe Marchi. Paresce. fc Pulonei (|2007l . hereafter DMPP) 
reported a surprising relation (see Fig. [TJ between the slope 
a of the low-mass stellar mass-function (dN/dm = m _Q! ) 
of globular clusters (GCs) and their concentration parame- 
ter c = log (rt/r c ), i.e. the logarithmic ratio of tidal- and 
core-radius. All high concentration clusters in their sam- 
ple had a steep mass-function (MF), while low concentra- 
tion clusters tended to have a natter MF (smaller a). The 
canonical IMF has a slope of a = +1.3 and a = +2.3 in 
the mass range 0.08 M - 0.5 M Q and 0.5 M - 150 M 



* e-mail: mmarks@astro.uni-bonn.de (MM); holger@astro.uni- 
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l|Kroupall200"ll : IWeidner fc Kroupall2004l ), respectiveljQ. The 
mechanism usually believed to be responsible for high cluster 
concentration is core-collapse driven by two-body relaxation 
and the same process causes low-mass stars to move to the 
outer cluster parts where they are removed by the external 
tidal field. So one would naively expect the relation to be 
the other way round, and this is indeed confirmed by the 
direct N-b ody simulations of init i ally n on-mass-segregated 
clusters of iBaumgardt fc Makinol (120031 ). 

Among other approaches to a solution discussed in 
DMPP, a variation of the IMF, e.g. depending on cluster 
density, could possibly solve this problem, if the present day 
concentration is primordial. This would for the first time be 
a strong observational evidence for a variable IMF. However, 
before adopting such an exciting solution, it is important to 
check for alternative solutions. 

Many observations show very young cluster s to be (par- 
tially) mass-segregated (e.g. recent pa pers by ISabbi et all 
l2007l : IChen et alj|2007l : iDib et al.ll2007l ). The reaction of ini- 
tially mass-segregated clusters to gas expulsion and the in- 
fluence which this has for the stellar MF in a cluster has 
never been investigated but may be important for under- 
standing the structure and stellar content of present day 
globular clusters. In this paper we study for the first time the 

1 Note: DMPP defined their MF as dN/dm = m+ a . 
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2 THE MODELS 



Figure 1. Observed trend between MF index a and the con- 
centration parameter c (taken from fig.l in DMPP). The dashed 
line is a fit to the distribution (Section 13.41 , The slope a is mea- 
sured in the mass-range 0.3 — 0.8 Mq and the canonical IMF has 
a rs 1.72 over this mass-range. 



influence of the gas expulsion process on the MFs of initially 
mass-segregated globular clusters (G Cs) with an invariable 
IMF using the recent results by iBaumgardt fc Kroupal 
|2007l . see Section^) for the dynamical response of star clus- 
ters to gas expulsion. 

The idea behind this is as follows: Star clusters form 
from massive, dense cores in giant molecular clouds. Stars 
form in these cores with a given star formation efficiency 
(SFE) reproducing the stellar IMF. The radiation of the 
newly born stars will immediately begin to disrupt their 
surrounding gas. Especially the winds of high-mass (O- and 
B-type) stars and the first supernovae will lead to ejection 
and disruption of the entire cloud. The stars have to react 
to the change in the potential and depending on the gas 
removal time-scale, the stars will or will not have time to 
react to this change, so by gas expulsion the final proper- 
ties of the cluster will be affected. The cluster expands and 
thereby preferentially stars from the cluster periphery are 
lost. If mass-segregation is a primordial feature of star clus- 
ters, mostly the low-mass stars are lost due to gas expulsion. 
Thus, the low-mass stellar MF changes strongly, while the 
high-mass part is left unaffected. Also the concentration is 
affected by cluster formation parameters such as the SFE 
and the gas expulsion time. Intuitively we expect the most- 
disruptive cases of gas expulsion to lead to remnant clusters 
that have the lowest density and largest depletion of low- 
mass stars. 

The paper is divided in three parts: In Section [2] we 
explain the models used and the procedure of our analysis. 
This is followed by a presentation of the results in Section [3] 
and a discussion & conclusion in Section [4] 



IBaumgardt fc Kroupal (|2007l . hereafter BK07) carried out 
a large set of TV-body integrations studying the effect of 
residual-gas expulsion o n the survival rate and final prop- 
erties of star clusters (cf. Lada. Margulis. fc Dearbornlll984 : 



iKroupa. Aarseth, fc Hurlevll200ll ). They placed star clusters 
on a circular orbit around a spherical galactic potential and 
varied (i) the star formation efficiency 



Mecl + Ma 



(1) 



where M GC i is the stellar mass of the embedded cluster and 
M gas is the mass of the residual gas within the cluster vol- 
ume, (ii) the ratio of the gas expulsion time to the crossing 
time of a cluster, TAf/£cross> and (Hi) the strength of the 
external tidal field in terms of the initial half-mass radius 
in units of the tidal-radius rjj/r*. In their models, they as- 
sumed that the SFE does not depend on the position inside 
the cluster, so gas and stars followed the same density distri- 
bution initially, which was given by a Plummer model. The 
influence of the gas on the stars was modelled as a modifi- 
cation to the equation of motion of stars. Gas expulsion was 
assumed to start at a c ertain time tp, which was set equal to 
one TV-body time unit (|Heggie fc Mathieulll986l ). equivalent 
to l/VE of a cr ossing time at the clusters virial (= gravi- 
tational) radius l|Binnev fc Tremainell 19871 ). After the delay 
time to, the gas density was decreased exponentially on a 
characteristic time-scale tm, the gas expulsion time-scale, so 
the t otal gas left at later times was given by l|Kroupa et al.l 
l200ll ) 



M gas (i) = M gas (0) e 



-(*-*d)/tm 



(2) 



All calculations were performe d with the collisional TV-body 
code NBODY4 (lAarsetblll999l ) on the GRAPE6 computers 
|Makino etaH 120031 ) at the Argelander Institute. All clus- 
ters contained 20000 equal-mass stars initially, distributed 
according to a Plummer sphere. On the one hand, the re- 
striction to equal-mass stars is necessary as otherwise a grid 
of cluster mass would need to be computed for each of the 
above parameter combinations. This is presently not feasi- 
ble. On the other hand, 20000 stars give good statistics. The 
integration proceeded for 1000 initial TV-body times (equiv- 
alent to about 300 initial crossing times). 

The above procedure captures the essence of the physics 
and their results can be summarised briefly as follows: Both 
the star formation efficiency and the speed with which the 
gas is removed have a strong influence on the evolution of 
star clusters. In the case of instantaneous gas removal (tm <C 
Across), clusters have to form wi th SFEs ^ 33 per cent in 
order to survive gas expulsion (|Boilv fc Kroupal l2003al lbT) . 
This limit is significantly lowered for gas removal on longer 
time-scales and clusters with SFEs as low as 10 per cent can 
survive gas expulsion in the adiabatic limit (tm ^> Across) 
if the external tidal field is weak. External tidal fields have 
a significant influence on the cluster evolution only if the 
ratio of r^/rt is larger than about 0.05. Below this value, 
star clusters behave nearly as if they are isolated. 
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2.1 Procedure & Assumptions 

For each of the 20000 stars in the models of BK07 we se- 
lected a mass from the can onical IMF using a C-routine 
descri bed in the appendix of iPflamm-Altenburg fc Kroupal 
(2006). This reproduces the steep MF of the high concen- 
tration clusters, which have values around a « 1.5 in the 
mass range 0.3 Mq — 0.8 Mq (hereafter referred to as the 
low-mass part of the MF) considered by DMPP, if MFs of 
these clusters still resemble their IMFfl We further omitted 
the brown dwarf (BD) part of the IMF, i.e. masses below 
0.08 Mq, since BDs have a negligible effect on the overall 
evolution of the cluster and they can't be detected in GCs. 

We distributed stellar masses among the 20000 stars 
according to their initial total energy: The star with the 
smallest specific energy, i.e. the star that is strongest bound, 
gets the largest mass and the one with the largest energy is 
assigned the smallest mass. However, a segregation of 100 
per cent as for this mass- assignment is most likeley exagger- 
ated and at least the outer regions of a cluster wouldn't be 
expected to be totally mass-segregated, so our assumption 
probes one extreme. We also performed our analysis with ini- 
tially not segregated clusters, where the masses have been 
assigned randomly to the stars to probe the other extreme. 
Note again, that the models of BK07 have equal-mass stars, 
but as the mass in gas in most cases is much larger than the 
total mass in stars, the dynamical evolution is dominated 
by the gas expulsion process for the time-span over which 
the models were computed. Since BK07 included the effect 
of an ext ernal tidal fiel d in the so-called 'near field approxi- 
mation' (|Aarsethll 19851 ). the specific energy of a star is given 
by 



E t , 



$(r) 



L 2 / Q 2 . 2 

-uj (6x e x + z e z 



(3) 



where the first term is the specific kinetic energy of a star, $ 
is the potential due to the residual-gas and the other stars 
in the cluster and the last term is a combination of cen- 
trifugal and tidal energy. The vector r = (x,y,z) is the 
position-vector of each individual star measured from the 
cluster centre, the e; (i — x, z) are unit vectors and the 
angular velocity of the cluster around the spherical galactic 
potential, 



GMg 



(4) 



is determined by the mass of the galaxy, Mg, and the galac- 
tocentric distance, Rg- 

Recently analytic techniques became avail- 
able to create m a ss-seg reg ation in modeled clusters 
llGaburov fc Gielesl 120081; ISubr. Kroupa. fc Baumgardd 



20081; iBaumgardt. de Marchi. fc Kroupal 1 20081 ). 

Gaburov fc Gielesl (12008) created a radius dependent 
'segregated MF' for the high-mass stars while the low-mass 
part is radius independent. While our procedure has all 
massive stars but no low-mass stars in the centre, their 
formulation of mass-segregation mixes both low and high 
mass stars in the inner part. Both techniques have no 
high-mass stars in the cluster outskirts. 



In order to determine the final concentration of the 
modelled clusters, the tidal- and core-radii are needed. In 
order to determine these, we first calculated the position of 
the cluster centre (projected centre in the case of th e core- 
radius) using the method bv lCasertano fc Hud |l985). 

The tidal radius of the cluster was then determined iter- 
atively, by first assuming that all stars still in the calculation 
are bound and calculating the tidal radius, 



n = 



GM CC 
3M G 



■Rg , 



(•>) 



at the end of the the iV-body runs, when the residual-gas is 
expelled. In a second step we computed the mass of all stars 
inside rt and used it to obtain a new estimate for the tidal 
radius from equation |(3J. This was repeated until a stable 
solution was found. 

To identify the core-radius, we placed radial annuli 
around the projected cluster-centre and searched for the ra- 
dius at which the projected number- density of stars dropped 
to half its central value. This method is of course sensi- 
tive to the choice of the width of the rings. A too small 
width leads to a very noisy surface-density profile, too large 
widths result in inaccurate values for the core-radius. Also 
the lower the number of stars at the end of the computations 
the noisier the surface-density-profile gets. So we decided to 
increase the ring-width for clusters with less than 5000 (25 
per cent) and 2000 (10 per cent) stars at the end of the 
A^-body integrations, respectively, in order to compensate 
for that effect. This hardly changes the measured core-radii 
for smooth density-profiles, but improves them significantly 
for those with a not so stable profile. The central-density 
was determined by averaging over the surface-densities in 
the first two to five annuli (depending on the ring-width) to 
reduce the influence of outliers. We took the mean value of 
the measured core-radii from three different projections as 
the best approximation of the true core-radius in the mod- 
els. Their standard-deviation was taken to be the error of 
r c . Models which survived gas expulsion but kept less than 
five per cent of their stars have been left out of consideration 
because radii-determinations were too difficult. 

With the tidal- and core-radius, the concentration and 
its error follow from 



c = log! 



ft 



(6) 



and error propagation. 

Furthermore, at the end of each integration, we 
plotted the MF of the stars that lie within the final 
tidal radius following the bin ning-method described by 
iMai'z Apellaniz fc Ubedal l|2005l . their experiment 3) assum- 
ing a power-law behaviour in the interesting mass-range: 
We chose to distrib ute the masses in 2 x N 2 ^ 5 bin s (second 
recommendation of iD'Agostino fc Stephens! 1 1 986 ) with an 
equal number of stars in each birjf]. Afterwards we deter- 
mined the slope a of the resulting MF in the mass range 
0.3 Mq — 0.8 Mq by applying a least-squares fit. Although 
the slope may change over this range (see upper right panel 
of Fig. [2}, we assigned just one a- value for each model. This 



2 A power-law fit to the canonical IMF gives a m 1.72 over this 
mass-range. 



3 Binning was performed over the whole mass-range 
0.08... 150 M . 
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Figure 2. Example mass-spectra for our analysis. The plots show the number of stars, N(m), per mass-interval normalised to the 
bin-width, Am, per bin versus the mass m, both on logarithmic scales. The left and right panels, respectively, show the same models (see 
title), one without (upper) and one with unresolved binaries (lower). The fraction of stars remaining in the cluster is / s t. Crosses indicate 
the IMF, open circles the resulting MF. The solid vertical lines limit the mass range 0.3 — 0.8 Mq. The dashed line is the fit to the final 
MF in this region. In mass-segregated clusters only the low mass-part of their MF is affected: The slope of the MF shifts to lower values 
for a if binaries are included in our analysis. Also note that the MFs with unresolved binaries are flatter for larger binary-fractions (cf. 
iKroupa. Gilmore. fc Toutlll99ll) . 



may not be completely representative for the actual mass- 
distribution for some models, but the uncertainties of the 
best-fitting line give an idea of how strong the data deviates 
from an ideal power law. The errors (indicated in Fig. [6} can 
be seen as the uncertainty in the determination of the slope. 



3 RESULTS 

3.1 The effect of gas expulsion on the stellar MF 

Fig. [3] depicts the effect of gas expulsion on the slope of 
the MF in the low-mass part in dependence of the cluster 
final concentration. As can be seen, the MFs of most high- 
concentration clusters still resemble their IMFs and don't 
show a strong change of the slope in the low-mass part, 
while for clusters with low-c several models do show a strong 
change in the slope. The exact position of a point in the dia- 
gram depends on the set of parameters for the corresponding 
model. 

In the upper part of the diagrams (above a « 1) one 



finds the models which kept most of their stars. These 
are preferentially models with a gas expulsion time-scale 
> icross- Except for the lowest concentration clusters, most 
of them still lie on their IMF slope. Those models with a 
smaller tm therefore have an extremely high SFE (above 
40 per cent) and experience a weak or intermediate tidal 
field only. Just three models with a gas expulsion time-scale 
larger or equal to the cluster's initial crossing time are found 
to be shallower than a w — 1 which is mainly due to their 
low SFE of < 15 per cent. Below a concentration of c ~ 0.8 
the strong tidal fields in these models lift even the large- tm 
clusters with a large SFE away from the IMF slope. 

In the left part of the diagrams (c < 1.4) several models 
show a strongly flattened MF (below a « 1), which are pre- 
dominantly the models with fast gas removal (tm <C t CIOSS ). 
These clusters experience intermediate or strong tidal-fields 
and have large SFEs (> 33 per cent), which is responsible 
for cluster survival. 

The upper right panel reveals the dependency of the fi- 
nal concentration on the the initial tidal field strength. Start- 
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Figure 3. The slope a of the MF at the end of the iV-body integrations in the mass-range 0.3 Mq — 0.8 Mq as a function of the 
cluster concentration c. The upper panels and the lower left show the same data, but coded for a different parameter used in the set of 
models of BK07: top left: SFE e; top right: strength of the tidal field measured by the ratio of r^/rf, bottom left: for the characteristic 



gas cxpulsion-time-scale, TMi measured in crossing-times, t c 



of the initial cluster. The horizontal dashed line gives the slope of the 



canonical IMF in the mass-range 0.3 Mq — 0.8 Mq; all initial models lie on this line and have 1 < c; n ; t < 2.5. The lower right plot 
compares the position of the points in the c-alpha-plane with (open symbols and crosses) and without (filled dots) binaries in dependence 
of the chosen binary-fraction. For the sake of clarity, error bars are omitted in these plots. Uncertainties are indicated in Fig. [6] 



ing almost exclusively with weak tidal fields at large con- 
centrations, stronger tidal fields become successively more 
common as the concentration decreases. This is easy to un- 
derstand: Due to general cluster expansion the tidal bound- 
ary of the cluster shrinks and the core radius grows, the 
final concentration can be expected to be lower for models 
in which the two radii are lying closer together initially, i.e. 
for the clusters in a stronger tidal field. Of course the fi- 
nal concentration is also affected by the gas expulsion time 



and the SFE, since they determine how large cluster expan- 
sion exactly is. That is, the models with different tidal field 
strength overlap slightly, but the overall c-value is estab- 
lished by the tidal field. Strong tidal fields at large galacto- 
centric distances may have occurred for those GCs forming 
in the immediate vicinity of pre-Milky Way gaseous building 
blocks. 

The lower left panel in Fig. [3] suggests a systematic 
behaviour of the MF slope with the concentration for a 
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Figure 4. Dependence of the slope a of the low-mass stellar MF (as in Fig. |3j on the gas expulsion time tm for a fixed value of the 
SFE e (as indicated in the plots). The horizontal dashed line shows the IMF slope. 



given gas expulsion time-scale. This can explicitly be seen 
for tm icross- The panels of Fig. [4] emphasize the strong 
dependence on this parameter for fixed SFEs. As the gas ex- 
pulsion time decreases the MF becomes shallower. Further- 
more, a stronger tidal-field leads to enhanced mass-loss and 
thus to a lower MF-slope. So the models with the strongest 
tidal-fields are the ones with the lowest value of a for a fixed 
gas expulsion time. The lower the SFE gets, the stronger is 
the effect of the gas expulsion time. For a high SFE, very 
strong deviations from the IMF slope start for the lowest gas 
expulsion times and the strongest tidal-fields first, migrating 
to larger tm and lower rh/r t as the SFE decreases. 

The influence of the exact set of parameters is so strong, 
that for gas expulsion times < 0.33 t CTOSS we even observe 
four strongly affected MFs in the high-concentration regime. 

3.2 Unresolved binaries 

An error occurring in any observational determina- 
tion of MF slopes is due to unresolved binaries. 
iKroupa. Gilmore. fc Tout] (|l99ll ) and lKroupal i|200ll ) showed 
that this leads to an additional apparent flattening of the 
MF. We tried to account for this fact in our analysis as fol- 
lows: We randomly picked stars from the stellar mass grid 
and randomly chose a secondary mass for them from the 
canonical IMF (random pairing). A 'centre-of-mass' star is 



then assigned the sum of the two selected masses. These 
'binaries' are still treated as single stars, i.e. they are 'un- 
resolved', but they now have the combined mass of their 
components. So the companions will stay together through- 
out the simulation and stay bound to or leave the cluster 
as a stellar system, i.e. possible hardening effects or disso- 
lution of binaries is not directly being considered. But it is 
included indirectly by assuming a variable binary-fraction 
in dependence of c. 

The binary-fraction of a cluster is defined as 

number of binary stars . 
m number of single stars + binary stars 

In order to account for a change of the binary-fraction of 
a cluster due to dynamical evolution, we vary /bin for dif- 
ferent models. In denser clusters close encounters between 
stars are more common than in less dense models (jKroupal 
Il995l ) . So dissolution of binary- and multiple-systems is more 
likeley in denser environments. Assuming that star clus- 
ters form with similar binary-fractions and that concentra- 
tion is a measure of cluster densitjQ, clusters with a high 
concentration after gas expulsion are expected to retain a 
smaller number of binaries than low-concentration clusters 

4 Depending on the actual value of rt and r c : in principle a cluster 
with larger c can be less dense than a cluster with lower c. 
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Figure 5. Left panel: The c — a plane for initially non- mass-segregated clusters with unresolved binaries. The dashed line is the relation 
proposed by DMPP (see Scction l3.4l l. Just a slight decrease of the slope can be seen due to enhanced binary-fractions in low-concentration 
clusters. Right panel: resulting mass-spectrum (as in Fig. [2} for no mass-segregation assuming unresolved binaries for the strongly affected 
cluster in the mass-segregated case with parameters e = 50%, tm = 0.33t C ross, = 0.15rt (compare with lower right panel of Fig. [J}- 
The shape of the MF stays more or less the same, i.e. all stars are equally affected by gas expulsion. 



after sufficient time for dynamical evolution. As has been 
shown in Section I3.ll the concentration at the end of our 
iV-body integrations approximately follows the strength of 
the tidal field. So we chose to distribute binary-fractions 
of /bin = 10%, 20%, . . . 60% according to models with 
tidal-field strength rn/rt = 0.01, 0.033, . . . 0.2, respectively. 
Thus, clusters experiencing a weaker tidal field are given a 
smaller binary-fraction and vice- versa. 

A choice of /bin independent from the concentration is 
not suitable to enhance the trend, because the strength of 
the shift depends on the binary-fraction. In principle one 
could explain the relation with unresolved binaries alone, if 
the true MFs of all clusters still resemble their IMF. But 
since GCs are old objects and should therefore be at least 
partially dynamically evolved, this is not expected for all 
of them (Section [1]). For our modelling, gas expulsion is in 
any case needed to initiate the trend and to span the whole 
range in concentration. 

After selecting the binaries and assigning masses we dis- 
tributed these stellar masses to the stars as described in 
Section [2. II and re-ran our analysis. 

As expected, unresolved binaries cause an additional 
shift of the MF slope. The effect can be seen in the mass- 
spectra of Fig. [2] (comparison of upper and lower panels) 
and its influence on the a — c-plane is shown in the lower 
right panel of Fig. [3] The slope of the MF changes accord- 
ing to the chosen binary-fraction towards lower values for a 
when compared with the mass-spectra or the location of the 
points in the a — c-plots, that consider gas expulsion only. 
The effect becomes stronger the larger the binary-fraction 
/bin is. By our assumption, the strongest change in the slope 
of the MF is seen in the least concentrated final models, be- 
cause these have been assigned the largest binary-fractions. 



All models with binaries are flatter than their IMF. The 
binary MF slopes differ from the non-binary MF slopes by 
Aa = |abin — a no -bin| = 0.005 (one of the /bin = 10 per 
cent models) to Aa = 2.8 (/bin = 60 per cent) with an 
average change of Aa = 0.5. 

3.3 Not-mass-segregated models 

Fig- (left panel) shows that primordial non-mass- 
segregated clusters do not show a strong trend of the MF 
slope with cluster concentration. Just a slight decrease of 
a is seen as c becomes smaller, which is due to the larger 
binary-fractions assigned for low c-values. 

The reason for this behaviour can nicely be seen in the 
mass-spectrum shown in the same figure (right panel). The 
diagram displays the same model parameters as Fig. [2] (lower 
right panel), but for the case of no mass-segregation. Here, 
the MF is just shifted along the y-axis while the slope hardly 
changes. Hence, stars of all stellar masses are equally re- 
moved from a non-mass-segregated cluster during expansion 
following gas expulsion leaving a unchanged. 

3.4 Comparison with observations 

The panels of Fig. [6] compare the influence of gas expulsion 
on the MFs of mass-segregated clusters with and without 
unresolved binaries, with the DMPP data. Fitting by eye, 
DMPP proposed a relation of the form (solid line) 

a(c) =-— + 2.5. (8) 

c 

In this Figure we added upper and lower limits (dashed lines) 
to their observational data 




Figure 6. Comparison of our results with the observational data by DMPP. The left panel shows the result considering gas expulsion 
only, the right one includes primordial binaries. The large open circles are the observed clusters, the small filled dots (with error bars) 
are from the iV-body integrations. The solid line is an eye-ball fit from DMPP (eq. ©); we added the dashed lines as limits to their 
observations (see text). The agreement is better for the models including the effect of unresolved binaries. 
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(c) = -— +2.5 

c 



a t (c) = -— +2.5. (9) 

c 



Both, the models with and without unresolved binaries, 
agree with the general trend that less concentrated clus- 
ters show a shallower MF, but the agreement is much better 
for the models including unresolved binaries. Even high-c 
clusters with a flattened MF are within or near the observa- 
tional limits. The not mass-segregated models, on the other 
hand, aren't able to reproduce the observed trend (Fig. [SJ). 
Unfortunately we do not reach values around c ~ 2.5 with 
our models, but we wouldn't expect significant deviations 
from the IMF there. Nevertheless, our data points suggest 
a slightly steeper run than that proposed by DMPP, which 
can be understood, because DMPP had reliable observa- 
tional data for only 20 GCs and the relation was an eye- 
ball fit. In addition, the N-body models of BK07 did not 
include the effect of dynamical evolution of the clusters be- 
yond the 300 initial crossing times (Section [3]), which will 
also ch ange the MF slope because of low-mass star evapo- 
ration (|Baumgardt fc Makinol |2003) . Dynamical evolution, 
however, also affects the concentration and the core-collapse 
process may produce concentrations as large as in the ob- 
served clusters with c ~ 2.5. The difference between the 
models and observed clusters may also be explained by the 
exact number of unresolved binaries. How clusters with dif- 
ferent initial stellar content and concentration evolve exactly 
after gas expulsion is still an open issue. 

The low-concentration model clusters in the upper 
left corner of Fig. [6] with only one observed cluster, 
both with and without binaries, are expected to shift 
to lower slopes during dynamical evolution. Such low-c 
clusters can indeed survive for a Hubble time and wan- 



der downwards in the c — a-plane to slightly larger con- 
centrations and join the strongly affected clusters be- 
fore they dissolve, if these cluster s start mass-segregated 
jBaumgardt, de Marchi. fc Kroupal 120081 , in prep.). The 
low-c clusters with a short gas expulsion time-scale, which 
expanded a lot and lost many low-mass stars, may be left in 
a collision-less state such that dynamical evolution doesn't 
affect the properties of these clusters much. 

Although we find several models (including binaries) 
with MF slopes —1 > a > —2 in the lowest concentra- 
tion regime, which follow the DMPP relation, they do not 
observe any clusters there. This may be understood be- 
cause most of the computed clusters in this domain expe- 
rience the strongest tidal fields and these models are un- 
likely to survive for a Hubble-time as they expand by an 
additional ~ 30 per cent due to stellar evolution and will 
probably dissolve rather qui ckly l|Fukushige fe Heggielll995l ; 
IBaumgardt fc Makinoll2003l ). 

The four model clusters at c ~ 1.7, which are more 
strongly depleted in low-mass stars than the other models 
in that part of the diagram are also further away from the 
observational data, although three of them are still within 
the observational limits set by us. This may be a selection 
effect, because we find just four out of 31 clusters above 
c = 1.4 with a strong effect on the MF slope and DMPP 
observed just six clusters in total in that regime. 

Note that we measured global MFs (GMF) from our 
integrations while in observations only local MF measure- 
ments are available, which are subsequently transformed to 
GMFs by assuming that the MF at the half-mass radius ap- 
proaches closely the GMF ( De Marchi. Paresce. fc Pulonel 



2000; IBaumgardt fc Makinol I2003T) . We checked this in our 
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N-body runs by looking at the local MF (LMF) of the stars 
near the half-mass radius, i.e. between the 40 and 60 per cent 
Lagrange radii, and found that we can reproduce the trend 
as well. The mean difference between LMFs and GMFs in 
our analysis is Aa = |o<gmf — chlmf! ~ 0.25. This difference 
might be removed by dynamical evolution. 



parameters may yield a similar MF slope (but also the con- 
centration is affected by them). This is especially true for 
clusters with less affected MFs, i.e. large tm and SFE, but 
also in other parts of the diagram some data points accu- 
mulate. In the case of ambiguous results, at least another 
parameter would be needed to lift this degeneracy. 



4 SUMMARY & CONCLUSIONS 

We have for the first time studied the effect of gas expulsion 
with and without unresolved binaries on the stellar MF of 
initially mass-segregated clusters with an invariable (canon- 
ical) IMF. 

Due to gas expulsion the cluster expands and stars in 
the cluster halo are more easily lost than stars near the 
centre. In a mass-segregated cluster this leads to the loss 
of the less massive stars causing a change of the low-mass 
MF. Our computations showed a flattening of the MFs in 
the low-mass part for less concentrated final clusters, but 
not for strongly concentrated clu s ters. This can explain 
the |Pe Marchi. Paresce. fc Pulonel (|2007T t data, while ini- 
tially non-mass-segregated models don't produce the ob- 
served trend. Depending on the parameters of the models, 
the final concentration and the slope of the MF are deter- 
mined. Primordial binaries induce an additional apparent 
shift in the MF slopes (compare panels in Fig. [6j| , which lead 
to an even better agreement with the DMPP data. If the dif- 
ference between theory and observation is due to unresolved 
binaries, adjusting the exact value of /bin of the model clus- 
ters to best reproduce the observational data could lead to 
predictions of binary- fractions for the observed clusters. 

The models can approximately be divided into two dif- 
ferent parts: The models with a high concentration of which 
most still have their IMF, and clusters with a low concen- 
tration, where a few systems lost a large portion of their 
(preferentially) low-mass stars and thus show a significant 
deviation from their IMF slope. 

Our models constitute a possible solution to the DMPP 
problem without invoking a variable IMF. This contribution 
emphasizes the importance of choosing physically realistic 
initial models for understanding cluster evolution: The abil- 
ity of our models to reproduce the observed c — a trend 
strongly suggests that initial mass-segregation and gas ex- 
pulsion play a critical role in star-cluster evolution and that 
we can't exclude the possibility of a universal IMF for all 
globular clusters. 

However, the MF slope and concentration in our 
calculations are measured shortly after the residual-gas 
has been expelled from the cluster once the surviv- 
ing clusters have re-virialized. Later dynamical evolution 
of initially mass-segregated clusters will additionally af- 
fect these two parameters. This topic is investigated in 
iBaumgardt, de Marchi. fc Kroupal (|2008l . in prep.). 

If gas expulsion is the determining factor for the ob- 
served trend between c and a then, by direct measurement of 
these parameters, observers may be able to constrain start- 
ing conditions of the clusters such as the SFE and the gas 
expulsion time-scale from Fig. [3] unless the c-a-plane is de- 
generate. From our results it is not clear if the slope of the 
MF is uniquely determined by just one set of parameters for 
a given concentration. A different combination of the three 
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